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PREFACE 

An  accurate  representation  of  the  physics  of  atmospheric  diffusion 
is  required  to  narrow  the  uncertainties  found  in  models  of  air  base 
defense  against  chemical  agents  delivered  by  tactical  ballistic 
missile.s.  This  Note  presents  a  discussion  of  the  diffusion  and 
dispersion  of  toxic  droplets  when  these  processes  are  affected  both  by 
atmospheric  turbulence  and  by  the  droplets'  own  considerable  fall 
velocity.  It  is  shown  how  currently  used  models,  which  may  neglect  the 
influence  of  fall  velocity,  can  easily  be  modified  to  include  this 
important  parameter. 

The  author  acknowledges  the  value  of  discussions  with  R.  Saucier, 

H.  Bach,  M.  Juncosa,  C.  Porter,  M.  Kamionski,  T.  Garber,  and  others.  He 
is  especially  grateful  to  S.  Hanna  and  W.  Krase  for  written  comments. 

This  work  was  performed  as  part  of  a  project  sponsored  by  the 
Defense  Advanced  Research  Projects  Agency. 
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SUMMARY 


This  Note  summarizes  a  preliminary  analysis  of  contamination 
footprints  associated  with  release  of  liquid  chemical  warfare  agents 
(C.W.)  in  the  atmosphere.  The  study  was  suggested  because  defense 
systems  analysts  often  require  contamination  footprints  for  situations 
that  depart  radically  from  the  limited  set  of  field  data.  Computer 
models,  presumably  based  on  a  sound  understanding  of  the  underlying 
physics,  are  then  used  to  predict  contamination  zones  and  magnitudes. 

We  found  that  the  C.W.  dispersal  and  diffusion  models,  developed  by 
the  Chemical  Systems  Laboratory  (CSL)  of  the  U.S.  Army  Armament  Research 
and  Development  Command  and  used  in  air  base  and  tactical  ballistic 
missile  defense  studies  by  Rand,  Institute  of  Defense  Analyses,  Hughes, 
General  Research  Corporation,  Aerospace  Medical  Research  Laboratory,  and 
others,  neglect  a  fundamental  element  in  describing  the  processes  of 
atmospheric  diffusion.  This  omission  translates  into  a  potentially 
severe  underestimate  of  C.W.  deposition  densities  and  a  physically 
implausible  picture  of  "toxic  cloud"  exposure.  Further  implications  and 
extensions  of  our  preliminary  findings  remain  to  be  explored. 

One  explanation  for  the  utility  of  the  CSL  computer  models  to 
defense  analysts  is  the  possibility  of  compensating  errors--that  the 
order  of  magnitude  error  in  modelling  diffusion  is  compensated  by  large 
omissions  or  cancelling  errors  in  the  other  components  of  the  models  in 
order  to  fit  a  limited  data  set.  If  this  is  the  case,  physical 
understanding  is  inadequate  but  the  results  may  be  satisfactory  for 
situations  that  are  close  to  those. of  the  few  field  tests.  Otherwise, 
the  models  fail  the  test  of  theoretical  plausibility  by  neglecting  the 
influence  of  fall  velocity  on  diffusion.  It  is  hoped  that  a  program  of 
theory  and  experiment  can  narrow  the  range  of  uncertainty. 
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I.  INTRODUCTION 


Defense  analysts  are  now  interested  in  understanding  the 
implications  of  a  potential  tactical  missile  threat  that  includes 
warheads  containing  chemical  warfare  (C.W.)  agents  intended  for 
atmospheric  dispersal.  Quantitative  understanding  of  this  threat,  for 
purposes  of  analysis  and  studying  defense  system  effectiveness,  relies 
heavily  on  our  ability  io  predict  the  contamination  tootpnnt  associated 
with  the  dispersal  of  such  agents.  Several  factors  that  influence  the 

biological  effects  of  these  agents  help  to  define  a  contamination 

footprint.  These  factors  include: 

•  Meteorological  conditions. 

•  Physical,  chemical,  and  rheological  properties  of  the  agents. 

•  Nature  of  the  agent  release  process  and  the  formation  of 

droplets  by  aerodynamic  stripping  as  influenced  by  dispense 
altitude  and  velocity. 

•  Formation  of  vapor  clouds,  and  the  tradeoff  between  liquid 
deposition  and  evaporation. 

•  Atmospheric  dispersal  and  diffusion  of  droplets  and  vapor. 

Because  field  data  and  observations  are  sparse,  computer  models  of 
the  underlying  physical  processes  must  be  used  to  develop  a  picture  of 
the  dispersion,  diffusion,  and  surface  evaporation  phenomena  that 
influence  the  contamination  footprint. 

Accurate  models  of  the  atmospheric  diffusion  of  toxic  clouds  of 
liquid  C.W.  droplets  should  be  key  ingredients  in  programs  to  compute 
contamination  patterns  and  C.W.  deposition  rates.  Such  models  have  been 
developed  by  the  British  [1]  over  the  past  twenty-five  years  in  the 
context  of  research  on  chemical  defense.1  Unfortunately,  the  widely 
used  U.S.  model  NUSSE2  [2],  developed  by  the  Chemical  Systems  Laboratory 
(CSL)  of  the  U.S.  Army  Armament  Research  and  Development  Command  fails 

'See  the  appendix  for  a  review  of  the  limited  theoretical  and 
empirical  data  base. 
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to  incorporate  the  British  results  in  its  current  version,  and  as  a 
consequence  overestimates  atmospheric  diffusion.  This  overestimate 
yields  C.W.  deposition  rates  and  contamination  footprints  that  could  be 
in  error  by  an  order  of  magnitude  or  more.  As  discussed  below,  we  have 
made  a  p:ri,,iinary  study  of  the  potential  discrepancy.  Figure  1 
exhibits  the  magnitude  of  the  possible  error  in  terms  of  the  deposition 
of  an  agent  released  from  a  point  source  at  a  height  H  above  the  ground. 
The  lower  curve,  A,  is  a  nondimens iona 1  representation  of  the  deposition 
rate  based  or:  an  atmospheric  diffusion  theory  that  is  similar  to  chat 
used  in  NUSSE2  and  other  CSL  programs.  The  upper  curve,  B,  is  based  on 
an  atmospheric  diffusion  theory  that  incorporates  a  more  plausible 
representation  of  the  physics  of  atmospheric  spreading  of  clusters  of 
liquid  droplets.  There  is  a  large  difference  in  the  deposition  pattern: 
the  theoretically  correct  pattern  involves  much  higher  rates  of 
deposition  (over  smaller  crosswind  distances).  The  difference  between 
the  two  curves  stems  from  the  following: 

Curve  A,  the  lower  one,  assumes  that  liquid  droplets,  as  they  fall 
through  the  atmosphere,  randomly  diffuse  like  parcels  of  air  about  their 
centroid--their  diffusion  characteristics  are  assumed  to  be  the  same  as 
tracers  with  air-like  mass,  and  are  entirely  influenced  by  the  turbulent 
eddy  motion  of  air  parcels. 

Curve  B,  the  upper  one,  reflects  the  finite  m3ss  of  the  falling 
droplets  to  the  extent  that  their  considerable  terminal  fall  velocity 
inhibits  the  random  diffusive  motion  about  their  centroid.  Thus  C.W. 
droplets,  as  they  fall,  are  only  slightly  buffeted  by  the  air's  eddy 
motion  and  in  fact  tend  to  fall  through  or  cross  the  trajectory  of 
eddies  rather  than  follow  them. 

C.W.  droplets  that  are  large  enough  to  evaporate  slowly  and  to 
survive  until  ground  impact  tend  to  be  between  one-half  millimeter  and 
five  millimeters  in  diameter,  and  have  terminal  vertical  or  fall 
velocities  between  1  and  10  m/sec.  For  these  droplets,  the  true 
diffusion  process  is  up  to  an  order  of  magnitude  weaker  than  one 
computed  by  neglecting  the  fall  velocity  effect. 

It  is  possible  that  the  CSL  models  employ  diffusion  oarameters  that 
are  thus  unrealistically  high  in  order  to  account  for  droplet 
depositions  actually  oDserved  during  the  limited  field  tests.  This 
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Fig.  1  -  Deposition  along  centerline 

implies  that  other  phenomena,  either  omitted  or  not  accurately  described 
by  the  CSL  models,  were  present  during  the  field  tests,  and  computer 
models  rely  on  the  "theory  of  compensating  errors."  Particle  inertia, 
wind  shear  and  the  turning  of  the  wind  vector  with  altitude  (Ekman 
spiral)  are  examples  of  phenomena  which  are  known  to  influence 
diffusion,  and  which  are  not  accounted  for  in  the  CSL  model.  However, 
none  of  these  mechanisms  are  as  important  as  the  vertical  velocity 
effect  over  the  time  scales  associated  with  hypothetical  C.W.  release 
altitudes.  Furthermore,  all  diffusion  processes  are  expected  to  be 
strongly  damped  by  the  fall  velocity  effect. 
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II.  DISPERSION  AND  DIFFUSION 

Atmospheric  winds  and  turbulence  cause  clouds  of  liquid  droplets  to 
disperse  and  diffuse  in  two  different  ways. 

1.  Because  the  clouds  include  droplets  of  varying  size  and  fall 
velocity,  they  tend  to  disperse  ballistically  in  the  downwind 
direction  as  drops  fall  and  are  carried  by  the  wind.  Small 
droplets  travel  further  downwind  than  do  larger  droplets.  This 
can  be  computed  from  droplet  size  distribution,  wind  speed,  and 
direction.  This  is  a  deterministic  process  that  we  designate 
as  ballistic  dispersion,  to  distinguish  it  from  diffusion. 

2.  The  crosswind  spreading  is  entirely  due  to  the  diffusion 
effect--even  if  all  droplets  were  a  single  size  they  would 
still  exhibit  random  motion  and  diffuse  crosswind  due  to 
atmospheric  turbulence.  They  would  also  diffuse  isotropically 
about  their  ballistic  trajectories.  This  can  be  accounted  for 
by  considering  a  polydisperse  cloud  of  particles  consisting  of 
subclouds  of  uniform  particles,  each  with  a  distinct  fall 
velocity.  The  entire  cloud  is  then  accounted  for  by  summing 
over  the  contribution  due  to  the  individual  subclouds.  Thus 
the  diffusion  process  is  linked  to  the  dispersion  process,  and 
both  must  be  described  accurately. 

If  we  consider  a  subcloud  (or  puff)  of  monodisperse  droplets 
released  at  time  t  =  0,  the  concentration  of  agent  is  represented  by  a 
Gaussian  spatial  distribution  with  a  single  standard  deviation  o(t)  in 
all  three  directions,  assuming  isotropic  turbulence  and  that  the 
"falling  droplet"  effect  is  felt  in  x,  y,  z  directions,  or 


(x  -  xQ(t))2  +  (z  -  zo(t))2  +  y2 
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where  N  is  mass  released,  x  is  distance  downwind  measured  from  release, 

z  is  vertical  distance,  y  is  crosswind  distance,  and  x  (t),  z  (t),  0  an 

o  o 

the  coordinates  of  the  centroid  of  the  cloud.  Note  also  that 


ft 

x  ( t 'I  =  /  udt  where  u  is  the  horizontal  windspeed, 

°  J  o  ft 

z  (t)  =  H  -  /  vdt  where  v  is  fall  velocity, 

°  Jo 

and  H  is  the  height  of  the  release. 
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Smith  and  Hay  (31  developed  a  theory  that  accounts  for  the 
dependence  of  o  on  finite  fall  velocity.  (See  appendix  for  further 
discussion.)  The  theory  is  surprisingly  tractable,  but  does  not  yield  a 
simple  explicit  formula  for  o(t).  Thomas  [4]  suggested  that  the 
Smith-Hay  theory  for  diffusing  falling  particles  can  be  simplified  to 
derive  the  simple  approximation  (for  preliminary  applications  and 
studies ) : 


•  • 


(2/3) 


V  v2  +  u2/62 


(2) 


where  o 
i 
t 
v 
u 
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is  the  standard  deviation  of  the  cloud 

is  the  turbulence  intensity  (magnitude  about  .1) 

is  the  time  since  release 

is  the  fall  velocity 

is  the  wind  speed 

is  the  ratio  of  Lagrangian  and  Eulerian  scales 
(magnitude  about  5). 
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III.  PRELIMINARY  ANALYSIS 


The  purpose  of  this  preliminary  analysis  is  to  explore  the  effect 
on  the  predicted  deposition  rate  when  the  large  fall  velocities  of  C.W. 
droplets  are  included  in  the  diffusion  process.  Therefore,  we 
concentrate  entirely  on  the  dispersal  and  diffusion  process  to  develop 
analytic  formulas  for  purposes  of  comparison.  Phenomena  such  as  droplet 
evaporation,  finite  initial  cloud  dimensions,  altitude  variations  in 
fall  velocity  and  wind  speed,  vapor  formation,  surface  evaporation,  and 
dosage-time  computations,  which  are  included  within  the  NUSSE2  model, 
are  best  analyzed  via  modifications  to  existing  computer  programs.  It 
is  not  difficult  to  modify  the  CSL  program  to  reflect  a  truer  picture  of 
the  diffusion  process.  That  portion  of  the  NUSSE2  program  dealing  with 
surface  evaporation  may  also  require  modification.  As  suggested 
earlier,  an  unrealistically  large  diffusion  parameter  could  (because  of 
other  errors  or  omissions  in  the  model)  be  obtained  as  a  result  of 
fitting  NUSSE2  results  to  limited  experimental  data,  but  we  have  rot  yet 
analyzed  the  limited  field  data  in  sufficient  detail  to  test  this 
hypothesis . 

THE  GAUSSIAN  PANCAKE  MODEL 

When  a  cloud  is  released,  it  contains  a  polydisperse  mixture  of 
droplets  of  different  diameters  and  hence  different  vertical  velocities. 
For  convenience,  we  use  a  distribution  of  terminal  fall  velocities  N(v), 
where  MQN(v)dv  is  the  mass  contained  in  the  velocity  interval,  v,  v  + 
dv.  A  cloud  thus  consists  of  a  continuum  of  subclouds,  each  with 
distinct  fall  velocity  v.  For  concreteness,  we  choose  a  log-normal 
distribution  for  N(v),  with  specified  median  velocity  and  standard 
deviation,  but  this  is  not  essential.  (Another  approach,  indicated 
below,  is  to  employ  a  distribution  for  drop  size  and  then  relate 
terminal  fall  velocity  to  drop  size.) 

For  a  polydisperse  Gaussian  cloud  or  puff  released  at  time  t  =  0  at 
x  =  0,  y  =  0,  »  =  H,  the  concentration  is 
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If  the  wind  speed  is  constant,  xQ(t)  =  ut ,  and  if  the  specific  fall 
velocity  is  v,  then  ®&(t)  =  H  -  vt .  (Note  that  a  depends  on  both  v  and 
t.) 

Let  D(x,y)  —  total  amount  of  liquid  deposited  per  unit  area  at  a 
point  on  the  ground  x,  y.  If  vt  >  o,  or  H  >  o,  then  the  surface 
concentration  does  not  change  appreciably  as  the  cloud  deposits  on  the 
surface.  Physically  it  may  then  be  shown  that  the  total  amount 
deposited  on  the  element  of  surface  area  dx  •  dy  is  equal  to  the  liquid 
contained  in  a  volume  element  whose  base  is  dx  •  dy  and  which  extends 
from  the  bottom  of  the  cloud  to  the  top.  Thus,  each  subcloud  with  a 
specific  velocity  v  deposits  liquid  droplets  as  if  it  were  a  thin 
Gaussian  pancake,  with  a  deposition  per  unit  area  given  by 


M  •  N(v)dv 
o 

2no2 


exp 


where  x  (v)  is  the  x-location  of  the  centroid  of  the  subcloud  of  fall 
o 

velocity  v  as  it  impacts  the  ground.  If  the  cloud  is  released  at  height 

H,  time  t  =  0,  and  distance  x  =  0,  then  the  subcloud  impacts  the  ground 

at  time  t  =  H/v,  and  the  x  location  of  the  this  impact  is  xc(v)  =  u*H/v. 

We  choose  u  as  a  characteristic  wind  speed,  taken  as  1/H  /"  udy,  the 

Jo 

average  speed  between  0  and  H.  This  Gaussian  pancake  model  is  analogous 
to  a  tilted  plume  and  to  the  model  employed  in  Ref.  2. 

As  discussed  earlier,  the  standard  deviation,  o,  is  approximately 
related  to  v  and  t  through 
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and  if  t  is  replaced  by  H/v,  we  obtain  the  specific  o(v)  for  the  sub¬ 
cloud  of  velocity  v. 


o(v)  =  3- 


i2  •  u2  »  H/v 

VTTTw- 
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SPATIAL  DEPOSITION 

The  spatial  deposition  D(x,y)  is  obtained  by  integrating  over  the 
depositions  from  each  subcloud  or 


D<x*y>  =  2^ 
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N(v) dv 
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If  o/H  <  1,  this  may  be  approximated  by1 


2  u  •  N/—  \  •  exp 

1-  W*)2I 

D(5.v)“-  x  . 

'o  > 

V  2n  •  x  *o 


(7) 


where  D(x,y)  •  H  /Mo  is  the  nondimensiona 1  deposition  rate,  x  =  x/H, 
y  =  y/H,  o  =  o/H,  and 


o  = 


-  (2/3)  •  iz  •  x  -  6 


^1  +  02/x2 


(8) 


For  a  log-normal  distribution  in  fall  velocity,  N(v)  or  N(u/x)  may  be 
written 
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‘This  formula  is  identical  to  the  Porton  model  for  the  deposition 
of  large  droplets,  except  for  our  specification  of  o  in  terms  of  fall 
velocity . 
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N(u/x) 


(u/x) 


exp  |-  [i.n(u/x)  -  y]2/2v2j  (9) 


where  exp  (y)  is  the  median  fall  velocity  and  v  is  the  standard 
deviation  of  fcn  v. 

To  calculate  the  deposition  rate,  several  parameters  must  be 
specif ied--i~wo  for  the  distribution  of  fall  velocities  y  and  v,  and 
three  for  the  atmospheric  diffusion  process,  u,  B  and  i.  Pasquill  (as 
well  as  others)  suggests  that  3  (the  ratio  of  the  Eulerian  to  Lagrangian 
lengths  scale)  is  between  1  and  10,  with  5  being  a  widely  quoted  value. 
He  also  suggests  that  a  turbulent  intensity  level  i  of  10  percent  is 
reasonable  for  most  atmospheric  conditions.  For  preliminary 
comparisons,  it  is  then  acceptable  to  set  3  at  5  and  i  at  .1.  Tne 
current  edition  of  Pasquill's  classic  treatise  (1)  suggests  that  &i  may 
be  roughly  constant  with  a  value  close  to  .6.  Presumably,  f5  depends  on 
meteorological  conditions. 

If  v  is  terminal  fall  velocity  arid  an  empirical  function  6(v)  is 
used  to  specify  droplet  diameter,  6,  in  terms  of  fall  velocity,  then  the 
mass  distribution  in  velocity  space  may  be  rewritten  in  terms  of  droplet 
diameter  or 


N(v)  =  Q(S)  *  ~ 


when  Q ( 6 )  is  the  mass  distribution  function  in  droplet  diameter  space. 
Given  representations  of  Q ( 6 )  and  6(v),  it  is  easy  to  determine  the 
deposition.  To  illustrate,  when  Q(6)  is  log-normal  in  droplet  diameter 
6  with  the  standard  deviation  of  £n6  equal  to  a,  and  the  median  droplet 
diameter  equal  to  exp(b) ,  the  deposition  relation  now  becomes 


D(X  y)  .  «i  =  i_  .  i  .  i-  .  exp|_l(4+  (£.nli)-b22 

V  ,y;  M  211  x  oa  dUnv  P  2  -2  2 

o  V  a  a 


.nWXV.V.  . 


*>  v  •_* 


*-*  *-•  v  *.*  *-?  v  v ^  r j  v 


11 


Since  the  downwind  motion  of  particles  is  assumed  to  be  ballistic,  the 
relation  v  =  u/x  translates  into  a  unique  relation  between  6  and  x,  or 
6(v)  =  6(u/x).  Thus,  the  computation  of  D(x,y)  requires  6(v),  the 
magnitudes  of  a  and  b  to  define  the  log-normal  distribution,  and  Eq .  (8) 
which  defines  o. 

The  illustrative  res  ts  shown  in  Figs.  1  and  2  are  based  on  Eq . 

(9)  rather  than  Eq.  (10).  The  curves  are  thus  suitable  for  comparing 
relative  deposition  densities  predicted  by  either  neglecting  or 
including  fall  velocity  Estimating  the  correct  absolute  magnitude  of 
the  deposition  density  would  presumably  require  6(v)  end  the  use  of  Eq. 
(11).  As  described  below.  Figs.  3  and  4  account  for  this  in  an  ad  hoc 
way  • 
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IV.  RESULTS  AND  DISCUSSION 


GENERAL 


We  have  computed  deposition  patterns  for  a  variety  of  wind  speeds 
and  distributions  of  fall  velocity.  Equations  (7),  (8),  and  (9)  form 
the  bases  for  Fig.  1 ,  which  indicates  the  deposition  along  the 
centerline  (y  =  0)  when  u  =  5  m/sec;  the  terminal  velocity  is 
characterised  by  a  media;;  of  5  m/sec  with  a  standard  deviation  of  3.5 
m/sec.  As  observed  earlier,  the  discrepancy  between  curves  A  and  B  is 
large  (over  an  order  of  magnitude).  Roughly  speaking,  and  when  all 
other  parameters  are  constant,  the  ratio  of  the  "true"  maximum 
deposition  rate  for  heavy  droplets  to  one  computed  using  the 
inappropriate  tracer  particle  diffusion  value  is  (along  the  downwind 


direction  y  =  0) 


t  27-2 

1  +  6  /x* 


Thus,  the  true  deposition  rate  is  always  greater  than  one  computed 
using  the  "light"  particle  diffusion  coefficient.  Conversely,  mass 
conservation  requires  that  the  true  crosswind  dimension  of  the 
contaminated  region  be  correspondingly  less  than  one  estimated  using  the 
incorrect  diffusion  theory.  This  is  demonstrated  in  Fig.  2,  where 
several  contours  of  constant  spatial  deposition  density  are  shown. 
Contour  A  corresponds  to  a  nondimensional  deposition  of  five  units 
computed  using  the  incorrect  NUSSE2  theory,  and  Contour  B  corresponds  to 
the  same  rate  of  five  units  when  computed  using  the  more  accurate 
theory.  Note  the  difference  in  crosswise  dimension.  .he  "true" 
contamination  region  is  narrower  than  the  "in'  rrect’  region. 

An  even  more  interesting  and  disconcerting  discrepancy  arises  when 
we  consider  a  contour  corresponding  to  a  deposition  of  20  units.  Such  a 
contour  (Contour  C)  may  be  found  using  the  correct  diffusion  theory,  but 
predicted  deposition  rates  are  always  less  than  this  value  when  the 
incorrect  theory  is  used.  Thus,  the  diffusion  theory  used  in  the  CSL 
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computer  program  always  underestimates  the  magnitude  of  liquid 
deposition,  and  yields  low  values  for  maximum  exposure  to  toxic 
droplets.  The  practical  importance  of  this  is  given  by  a  hypothetical 
example:  suppose  the  nondimensional  value  of  20  units  corresponds  to  a 

biological  threshold  effect  for  a  C.W.  release  of  500  kg  at  10  km.  The 
incorrect  theory  then  results  in  a  significant  underestimate  of  hazard 
to  personnel  from  liquid  contamination.  However,  as  we  have  speculated, 
although  individual  elements  of  C.W.  defense  models  may  be  in  serious 
error,  the  overall  output  cf  the  models  might  correspond  to  test  data 
obtained  in  limited  field  tests.  This  point  deserves  further  study. 


APPLICATION  TO  POWER  LAW  MODEL 

For  the  sake  of  simplicity  and  theoretical  consistency  we  have  used 
a  relation  for  o  based  on  Smith  and  Hay's  theory  for  the  diffusion  of  a 
falling  monodisperse  cloud  or  puff.  (See  appendix  for  further  details.) 
Because  there  are  few  data  about  such  clouds,  and  because  turbulent 
diffusion  theory  is  limited  to  idealized  situations,  some  investigators 
prefer  to  use  empirical  diffusion  coefficients  and  o's  similar  to  those 
that  describe  the  spreading  of  continuous  plumes,  such  as  smoke  or 
particles  issuing  continuously  from  an  elevated  chimney.  (As  described 
in  the  appendix,  the  plume  diffusion  process  differs  from  cloud  or  puff 
diffusion  except  at  very  long  times.)  In  this  spirit,  NUSSE2  uses  the 
plume  formula  o(x)  =  A  •  x*,  where  A  and  a  are  empirical  coefficients, 
presumably  derived  from  the  spreading  of  plumes  of  iight  tracer 
particles.  In  the  spirit  of  Thomas'  simplification  of  the  Smith-Hay 
results  (4),  we  suggest  that  a  fall  velocity  correction  for  a  puff  would 
lead  to 


o(x,v)  =  A 


This  may  be  further  simplified,  when  o/H  <  1 ,  to 


v  '.w.v  v.V. v  Wc  -  o ,'i ,  .v. 
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where  0  is  again  assumed  to  be  close  to  5.  This  suggests  a  simple  way 
to  crudely  correct  for  vertical  velocity  in  the  NUSSE2  results.  A  more 
rigorous  correction  procedure  is  only  slightly  moie  complicated. 

Figures  3  and  4  demonstrate,  in  a  practical  context,  the  magnitude 
of  the  discrepancies  between  the  correct  and  incorrect  diffusion 
theories.  Consider  a  release  of  500  kg  at  1.50  km.  If  the  average  wind 
speed  is  7.1  m/sec,  the  downwind  distribution  of  liquid  deposition  is 

maximum  at  about  1.1  km.  However,  the  corrected  maximum  deposition  rate 

4  2  3 

is  1.5  10  mg/mm  ,  whereas  the  uncorrected  maximum  is  2.5  x  10  --roughly 

a  factor  of  5.6  (see  Fig.  3).  The  crosswind  effect  is  shown  in  Fig. 

4--the  crosswind  distance  of  appreciable  deposition  shrinks  from  about 

700  meters  to  120  meters ‘--the  deposition  region  shrinks  and  closely 

simulates  a  line  source  of  material  for  subsequent  evaporation. 


1.5  x  104 

1.0  x  104 

0.5  x  10“ 

0.5  1.0  l.b  2.0  2.5  3.0  3.5 

Downwind  distance,  x,  km 

Fig.  3  -  Centerline  deposition  based  on  ad  hoc  correction 
‘From  symmetry,  total  crosswind  distance  =  2y. 
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indicated  by  Pasquill  and  Smith,  plume  diffusion  is  expected  to  be  more 
weakly  affected  than  puff  diffusion,  and  the  limited  effects  that  are 
observed  are  quite  consistent  with  theory.  Furthermore,  the 
microspheres  used  in  the  Canadian  tests  are  considerably  smaller  than 
the  probable  droplet  sizes  that  would  result  from  the  atmospheric 
dispersal  of  liquids  carried  by  tactical  missiles. 

The  discussion  and  analysis  described  in  this  Note  are  based  on  a 
decades  old  theory  that  has  never  been  adequately  tested  against  data. 
Because  there  is  now  a  revival  of  interest  in  C.W.  defense,  ar.d  because 
the  puff  diffusion  process  for  large  droplets  differs  considerably  from 
the  diffusing  of  plumes  of  aerosols  or  atmospheric  pollutants  important 
in  environmental  studies  (which  have  been  studied  intensively  in  recent 
years),  we  believe  that  now  is  the  time  to  develop  an  appropriate 
experimental  and  analytic  program.  Such  a  program  would  include 
"research  grade"  experiments  (using  microspheres,  or  similar  particles 
of  known  properties)  and  field  tests,  performed  carefully  and  based  on 
Theologically  appropriate  simulants. 

A  suitable  program  of  test  and  analysis  would  provide: 

•  Comparison  with  experiment  and  field  data. 

•  Confirmation  of  the  diffusion  theory  discussed  here  for 
appropriate  fall  velocities  using  both  alternative  theories  and 
Monte  Carlo  simulation.  Dr.  Steven  R.  Hanna  has  suggested 
(personal  communication)  that  a  Monte  Carlo  diffusion  exercise 
could  be  conducted  using  velocity  fields  generated  by  a  Large 
Eddy  Simulation  model,  currently  under  study  by  an  Army 
Research  Office  working  group. 

•  Consideration  of  the  impact  of  droplet  inertia  on  atmospheric 
diffusion. 

•  Consideration  of  the  liquid-water  evaporation  process  on  the 
ground,  and  the  vapor  exposure  problem. 

•  Consideration  of  meteorological  effects  which  are  particularly 
important  when  releases  within  and  above  the  atmospheric 
boundary  layer  are  analyzed.  (The  CSL  model  employs  a  wind 
profile  which  is  inappropriate  for  release  above  several 
hundred  meters.)  The  influence  of  wind  shear  and  the  variation 
of  wind  direction  with  altitude  should  be  included. 
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We  remark,  in  closing,  that  such  a  program  would  both  clarify  the 
physics  of  droplet  diffusion  and  narrow  the  uncertainty  in  models  of 
C.W.  defense. 
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Appendix 

SURVEY  OF  PAST  WORK  ON  THE  DIFFUSION  OF  FALLING  CLOUDS 

PUFF  AND  PLUME  DIFFUSION 

Studies  of  atmospheric  diffusion  distinguish  between  plume 
diffusion  and  puff  diffusion.  Plume  diffusion  refers  to  circumstances 
where  material  release  and  sampling  times  arc  long  compared  with  travel 
time  from  the  source,  and  puff  diffusion  refers  to  circumstances  where 
material  release  and  sampling  times  are  short  compared  with  travel  time. 
Thus,  an  instantaneous  source  leads  to  puff  diffusion  and  a  continuous 
source  to  plume  diffusion.  The  transport  and  diffusion  of  a  cloud  of 
material  suddenly  released  in  the  atmosphere  clearly  corresponds  to  a 
puff  diffusion  process.  The  physics  of  puff  diffusion  is  much  less 
precisely  characterized  than  that  of  plume  diffusion.  In  the  words  of  a 
recently  published  handbook  [5],  "We  have  a  few  theories  for  puff 
diffusion  and  a  data  set  that  is  several  orders  of  magnitude  smaller 
than  the  data  set  for  plume  diffusion."  Furthermore,  the  empirical 
literature  on  puff  diffusion  includes  only  limited  data  on  the  effect  of 
atmospheric  stability,  and  virtually  no  data  on  the  dispersion  of  clouds 
of  droplets  or  particles  with  appreciable  fall  velocity. 

Analysts  must  then  rely  on  a  blend  of  theory,  empiricism,  and 
judgment.  Despite  the  limited  data,  intuition  suggests  that  the 
diffusion  of  clouds  or  clusters  of  heavy  particles  must  be  inhibited  by 
the  "trajectory  crossing  effect."  This  effect  was  first  described  by 
Soviet  scientist  M.  I.  Yudinn  [6j:  "When  falling,  a  heavy  particle 
crosses  trajectories  of  air  particles  so  that  it  interacts  consecutively 
with  different  air  particles.  As  a  result,  the  succession  of  velocities 
of  a  heavy  particle  does  not.  coincide  with  individual  changes  of  tin 
velocity  of  an  air  particle."  This  leads  to  a  reduction  in  dispersion. 

Starting  with  Yudine,  a  number  of  investigators  have  analyzed  plume 
diffusion  of  heavy  particles,  but  only  Smith  [ 7  j  and  Smith  and  Hay  (3) 
have  written  about  the  puff  diffusion  counterpart.  Fuff  diffusion  is 
sensitive  to  a  narrower  range  of  turbulent  eddy  sizes  than  plume 
diffusion,  being  most  influenced  by  eddies  with  sizes  close  to  the  puff 
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size.  Similarity  laws  for  puff  diffusion  have  been  developed  by 

Batchelor,  but  such  laws  do  not  include  the  influence  of  fall  velocity. 

For  long  times  of  flight  of  light  particles,  the  Batchelor  relation 
2  3 

o  =  tt  ,  where  e  is  the  eddy  dissipation  rate,  has  been  recommended  for 
practical  use  (5]. 

The  puff  theory  for  heavy  particles  first  proposed  by  Smith,  and 
later  modified  by  Smith  and  Hay,  starts  with  Batchelor's  observation 
that  puff  diffusion  depends  on  the  relative  motion  of  pairs  of 
particles.  Tnus  a  two-particle  Lagrangian  covariance  is  required, 
whereas  plume  diffusion  relies  only  on  a  single-particle  Lagrangian 
covariance. 

Smith  and  Hay  assume  homogeneous  isotropic  turbulence,  invoke  the 
Taylor  transformation  x  a  Ut  to  relate  space  and  time  covariances,  and 
assume  the  validity  of  the  Hay-Pasquill  hypothesis  to  relate  the 
Lagrangian  time  covariance  to  the  Eulerian  time  covariance  in  terms  of 
B,  the  ratio  of  the  Lagrangian  time  scale  to  the  Eulerian  time  scale. 
Ultimately  the  Smith-Hay  equation  for  0,  in  terms  of  the  full  spectrum 
of  turbulence,  becomes 


do 

dt 


sinks 

ks 
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-o  k 


e. 
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dsdk 


(A.  1) 


where  E(k)  is  the  energy  density. 

If  Ut/B  is  large  enough,  then  further  analysis  (after  assuming  an 
exponential  correlogram)  leads  to 
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where  o’  =  o/t,  X1  =  X/2,  i  is  turbulence  intensity  and  l  is  the 
Eulerian  length  scale. 

If  Ut/B  is  not  large  enough  to  be  accurately  replaced  by  infinity, 
the  covariance  approach  leads  to 
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(A. 3) 


where  E  is  a  function  graphed  on  p.  92  of  Kef.  3.  Two  further 
approximations  for  o  are  also  developed: 

1.  If  the  initial  value  of  o'  is  sufficiently  small,  then  o'(0) 
may  be  set  equal  to  zero  during  the  integration  of  Eq.  (A. 2), 
and  a  universal  curve  is  obtained  for  initially  small  puffs  (or 
clusters)  .  This  curve  portrays  a  relation  y  =  F ( ) » 

where  X  is  measured  from  an  arbitrary  origin. 

2.  A  further  simplification  results  after  noting  that  E  in  Eq. 

(A.  3)  may  be  replaced  with  (reasonable  accuracy)  by  its  maximum 
value,  leading  to  Thomas '  approximation  (4). 
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(A. 4) 


LIGHT 

These  o's  correspond  to  o  ,  and  refer  to  the  diffusion  of  light 

particles . 


DIFFUSION  OF  HEAVY  PARTICLES 

The  results  discussed  above  were  originally  derived  for  clusters  of 
light  particles.  However,  they  may  be  extended  to  clusters  of  particles 
that  have  sufficient  mass  to  acquire  appreciable  terminal  velocity,  but 
are  still  small  enough  to  respond  to  the  local  turbulent  velocity 
fields . 

Under  these  conditions,  where  particle  inertia  may  be  neglected,  an 

approximate  equivalence  principle  is  described  by  Pasquill  and  Smith 

[1],  who  note  that  a  falling  particle  moves  horizontally  with  the  fluid 

and  experiences  Lagrangian  variations  equivalent  to  passing  through  the 

space  spectrum  with  velocity  U/3,  and  simultaneously  falls  through  the 

space  spectrum  with  velocity  v.  The  result,  they  suggest,  is  equal  to 

2  2  21/2 

moving  through  the  space  spectrum  with  velocity  (v  '  +  UVB  ) 


-VWV'V 


,  •  >  * 


22 


Frequencies  at  a  fixed  point  then  appear  to  be  increased  in  the 
2  2  21/2 

proportion  (v  +  U  /$  )  /U,  and  the  integral  time  scale  is 

correspondingly  decreased. 

When  v  is  nonzero,  the  parameter  6  in  Eqs .  (A.l),  (A. 2),  (A. 3),  and 

(A. 4)  may  be  replaced  by  &'  =  {(v/U)2  +  1/f)2}  .  For  example,  when 

o(0)  is  small,  o/l  may  be  represented  by 


Furthermore,  the  validity  of  replacing  6  by  &'  is  independent  of 
the  particular  form  of  the  correlation  function.  Thus,  Sawford's  recent 
study  of  relative  diffusion  (8],  using  alternative  representations  of 
the  two-particle  Lagrangiai:  correlation  function,  may  easily  be  modified 
to  include  the  influence  of  fall  velocity. 

The  heavy  particle  effect  is  greater  for  puff  diffusion  than  for 
plume  diffusion.  For  example,  Csanady's  plume  theory  [9]  leads  to 
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and  the  Smith-Hay  theory  (and  Thomas'  approximation)  for  puff  diffusion 
leads  to 
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For  droplets  which  fall  with  terminal  velocities  as  large  as  10 
m/sec,  we  see  that  the  standard  deviation  is  30  percent  of  the  light 
particle  value  for  a  plume,  and  10  percent  of  the  light  particle  value 
for  a  puff,  when  windspeed  is  5  m/sec  and  5  s  5. 

The  parameter  &  is  not  very  well  characterized--sorne  analysts  now 
believe  that  $i  is  a  constant  of  unit  order  (say  .6),  and  others  have 
tried  to  infer  a  dependence  of  6  on  stability. 

As  far  as  experimental  verification  is  concerned,  there  does  not 
appear  to  be  any  field  data  for  falling  clusters  that  may  be  used  to 
evaluate  or  modify  the  Smith-Hay  data.  Thomas  alluded  to  such  data  in 
his  1964  note,  but  we  have  not  yet  found  any  in  the  open  literature. 

Because  the  characteristic  dimension  o  for  falling  clusters  is  much 
less  than  for  clusters  of  light  or  tracer  particles,  we  speculate  that 
interactions  between  wind  shear,  either  horizontal  or  vertical,  and 
diffusion  are  less  important  than  for  light  particles. 

Particle  inertia,  and  its  effect  on  the  diffusion  of  falling 
droplet  clouds,  is  a  virtually  unexplored  area.  The  Smith-Hay  theory  is 
based  on  the  assumption  that  inertial  effects  are  not  important  except 
to  define  the  fall  velocity  of  the  particle.  Certainly,  we  expect  that 
inertia  is  negligible  for  droplets  in  the  submillimeter  range  (based  on 
plume  diffusion  results). 

We  are  less  certain  about  larger  droplets  in  the  1  to  10  mm  size 
range.  Analysis  by  Reeks  [10} ,  and  Nir  and  Pismen  (11),  indicate  that 
particle  inertia  has  little  effect  on  plume  (absolute)  diffusion  when 
the  turbulent  velocity  is  smaller  than  the  velocity  of  fall,  and  that 
Csanady's  formula,  (A. 5),  is  valid  over  a  broad  range  in  particle 
inertia . 

Thus,  although  particle  inertia  may  inhibit  the  large  effect  of 
fall  velocity  in  reducing  diffusion,  the  magnitude  appears  to  be  small 
for  plume  diffusion.  Until  the  influence  of  inertia  on  the  puff 
diffusion  process  has  been  analyzed  properly,  we  therefore  believe  that 
its  neglect  may  be  justified. 

In  summary,  then,  there  is  but  one  theory,  and  no  "research  grade" 
data  to  define  the  diffusion  of  clusters  of  falling  particles.  The 
theory  is  physically  plausible  and  indicates  an  important  diminution  of 
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dispersion  for  heavy  particles.  Verification  and  extensions  of  the 
theory  remain  to  be  pursued.  Clearly,  those  who  are  interested  in 
accurately  modelling  the  diffusion  of  toxic  droplet  clouds  should  now 
use  formulas  like  (A. 5)  and  (A. 6),  or  the  Pasqui 1 1 -Smith  described 
above,  to  modify  their  representations  for  a,  since,  as  shown  in  the 
body  of  the  Note,  the  heavy-particle  effect  can  lead  to  an  order  of 
magnitude  difference  in  deposition. 
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